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SECTION I 
INTRODUCTION 


This volume describes the operation of solidification Model 2 and It provides 
documentation of the software asisoclated with the model. Model 2 calculates 
the macrosegregation In a rectangular Ingot of a binary alloy as a result of 
unsteady horizontal axlsymmetrlc bidirectional solidification. A description 
of the physics and numerical techniques Is In Volume 1 of this report, 

The solidification program allows Interactive modification of calculation 
parameters as well as selection of graphical and tabular output. In batcii 
mode, parameter values are Input In card Image form and output consists of 
printed tables of solidification functions. The operational aspects of Model 
2 that differ substantially from Model 1 are described In Section 2. A complete 
guide for running Model 1 on the Prime ^00 system at MSFC Is In Volume 1 1 1 of 
reference 1. Section 2 assumes no familiarity with programming, but It does 
assume the reader is familiar with Volume 1. 

Section 3 of this report contains the global flow diagrams and data structures 
of Model 2. As with Model 1, the primary program documentation Is the code 
Itself. It is assumed that any reader requiring the Information In Section 3 
thoroughly familiar with the calculation desc-rlbed by Volume 1, with the program 
operating characteristics, and with the FORTRAN language. 


SECTION 2 
OPERATING GUIDE 


The operation of solidification Model 2 Is very similar to tha'. of Model 1. 

The input phase and output selection phase operate as for Model 1. The specific 
input parameters for Model 2 are described below. Output checkpoint control is 
a new procedure with Model 2; Interactive checkpoint control is described in 
Section 2.', and batch checkpoint control is described in Section 2.2. The 
mechanics of interacting with the Prime system are described In reference 1. 

2.1 RUNNING THE MODEL IN INTERACTIVE MODE 

The procedure for logging in to the Prime ^00 system, starting the model, and 
running the input and output phases is the same as for Mode* 1 with the exception 
of the commands for starting Model 2. After logging in, to invoke execution 
of Model 2 enter the commands 

CO MPS2 
SEG #MPS2 

Because of the time stepping nature of the unsteady model, the operation of the 
calculation phase Is unlike the operation of Model 1. After the input phase is 
complete, the model takes one step away from the chill and then prints the 
message shown below. 

CMTfR 

C 

This choice of action is given to the user at each checkpoint. Entering a G or 
a T wi 1 1 cause the program to go to the output phase in which the user can select 
specific graphical or tabular output from the calculation up to that point. 

When the user leaves the output phase, the selection shown above Is again 
presented. The user can switch freely between tabular and graphical output or 
ha can continue the calculation, start a new case or terminate execution. If 
the continue (C) option is selected, the program will display the following. 


T TO filSPLAV 
G TO BlSPLftV CRW»HS, 

C TO COHTINUE THIS CWE, 

0 TO TEKfllNATC RUN, OR 
N TO PROCEED TO NEXT CRSE. 
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TIK SINCE lEQINnING OF SOLIOIFICATIOH * I.470BIS ICC. 

LIOUIDUS ISOTHCRN IS •.7fl9995E>M CN FROM CHILL. 

EUTECTIC ISOTNERN IS *1.99299 CN FRON CHILL. 

ENTER EITHER THE TINE OF THE NEXT CHCCKFOIHT IN THE FORNRT T-UNLLf 

OR THE LIOUIDUS LOCRTION NT THE NEXT CHECICFOIHT IN THE FORNRT L'URLUl 
OR THE EUTECTIC LOCRTION RT THE NEXT CHECEFOIHT IN THE FORNRT I'URLUE 


The checkpoint procedure puts the user in control of the progress of the 
calculation. The stopping places specified by the user as checkpoints allow 
him to inspect the progress of the calculation or display the current state 
of the S/L zone at any point In the solidification process. For convenience 
the stopping point may be specified by time, by eutectic front location, or 
by liquidus isotherm location. For example, to set up the next checkpoint at 
the instant when eutectic forms at the chill, enter 

E«0. 

To stop the model when the dendrites first reach the centerline !n a mold of 
half-width 2 cm., enter 

L»2. 

To run the solidification through to completion enter a eutectic front position 
greater than or equal to the mold half-width. 

The unsteady model allows greater flexibility in plot setup than the stteady 
model with plot function scales brought under user control. In all cases the 
program will calculate a scaling and then give the user the alternative of 
proceeding with the calculated scaling or entering a different scale. For 
vector plots the scale is a single vector magnitude, for contour plots the 
variables are a base level value and an interval between levels, for profile 
plots the scale is specified as a range with a remote exponent and the number 
of tic intervals. If a nonzero remote exponent, say 3, is to be used the 

■3 

actual plot interval bounds should be multiplied by 10 before they are Input. 
Profile plots over the final solid also require that the user enter the x or 
y location of the profile, where x is always the distance from the chill in 
cm., and y is the distance from the botton of the ingot in cm. Any number of 
profiles can be put on one plot, but the line patterns will repeat after 6, 
Example ij.l shows plot scale input and profile selection. 




. 
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2.2 RUNNING THE MODEL IN BATCH MODE 

The batch operating procedure Is the same as for Model 1 wttt. the exception 
of some modifications to the Input files for specification of checkpoints. The 
command which starts execution of the model in batch mode is 

CO MP 52. BATCH 

Checkpoints are defined on cards immediately following the paramete" cards, 
one- checkpoint per card. The end of the case Is denoted by a card with the word 
END starting In column one. The checkpoint card has three fields. The first 
field begins in column 1 and must be one of the letters T, L or E denoting time, 
liquidus Isotherm position or eutectic front position. The second flfcld runs 
from column 6 to column 15 with the FORTRAN format E10.3» This field should 
contain the value specified for the checkpoint location. The. third field 
begins in column 16 and specifies the selection of tabular output to be printed 
at the checkpoint. Each column controls one table, with the assignment in the 
same order as the interactive tabular output selection shown in example A.l. 

A one In the column corresponding to a table will cause that table to be printed 
at the checkpoint, If the column Is blank, the table will not be printed. 

Thus, a one In column 17 would cause the final local average composition to be 
printed. Any number or combination of tables may be printed at eacli checkpoint. 

Batch cases can be stacked as before, with the case name card for case M 
following immediately after the END card of case M-1. A sample batch case is 
shown in Section k.2, 

2.3 INPUT PARAMETERS 

Described below are the parameters which define a particular case, of unsteady 
horizontal bidirectional solidification in a casting. The same set of parameters 
is used in both Interactive and batch modes. The default case Is the case 
that is built-in to the interactive I/O controller: The parameters take on 

their defaM'it values each time Interactive execution is initiated. There Is 
no default case for a batch run: all parameters must be specified for each 
case. The format for the batch cards Is given In Section 2.2 and reference 1. 
Interactive parameter selection and modification is described in reference 1. 
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2.3.1 Alloy Specification 

The binary alloy Is specified by chemical symbols exactly as Is done In Model 
1, and Model 2 uses the same data base as Model 1. The default alloy Is 
Al~^.5^Cu. 


3,3*2 Solidification Process Parameters 

The parameters in this group describe the conditions under which the casting 
Is made Including the mold geometry, the thermal conditions* and the strength 
of the gravitational force. 


Parameter Description 

Symbol 

Default 

LI mi ts 

Mold half width 

^c 

2.7 cm 

X >0 
c 

Mold height 

L 

6.35 cm 

L>0 

Isotherm time exponent 

q 

1 

q>0 

LIquIdus coefficient 


.017 cm/sec 

b^_>0 

Eutectic coefficient 


.017 cm/see 

b^>0 

Initial separation 


-2. cm 

a_<0 

t 

Gravitational force in 

g 

IG 

9>0 


-y dl recti on 


As noted in Section 3*2.2 of Volume I, the unsteady model is limited to cases 
in which the temperature of the chill Is brought below the eutectic temperature 
gradually rather than instantaneously. This restriction Implies a^ must be 
significantly less than the upper limit shown above. Results for cases with 
a^ close to zero are not predictable. 

2 . 3.3 Permeability 

The permeability model in the unsteady mcidel is identical to that In the steady 

""7 2 

state model. The defai.lt value of y is 6X10 ' cm . 


2 . 3 .^ Numerical Methods Control Parameters 


In the majority of cases the parameters in this group should 
from their default values, however operator control of these 
be useful In trouble-shooting. The parameters are explained 
of numerical methods In Volume 1 of this report and Volume I 


not be changed 
parameters could 
In the discussions 
of reference 1 . 
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Parameter Description 

Symbol 

Default 

Limits 

Number of horizontal 

^1 

20 

lHNj<50 

mesh points 
Number of vertical 


10 

3<Nj<50 

mesh points 


Maximum number of 

^SOR 

200 

«S0R>° 

pressure Iterations 


Pressure convergence 

‘^SOR 

10-^ 

E30,>5X10“ 

criterion 


Fraction of CFL step 

f 

.5 

0<f<l 

to take 




Maximum Step size 

AXg 

. 1 cm 

AXg>0 

In 



Initial at chill 

9U 

.9 

.Kg, I<1. 

Minimum horizontal mesh size 

'^imin 

A 

h<H. , <50 
— imln— 
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3.1.3 Input Section 

1 

CONTROLLER! 

INCON 


INITIALIZE INPUT ARRAYS 

CALLED FROM! 

INIT 

(INCON) 


INITIALIZE terminal 
CONTROL SYSTEM 
(lACTI) 


BATCH? 


READ CASE TITLE 
(BATCH I ) 


PUT ID PAGE ON SCREEN 
AND GET CASE TITLE 
(IACTI) 


^ yes / 

END OF FILE 7 > KSTOP 


GET ALLOY 
(IACTI) 


ACCESS ALLOY DATA BASE 
(ADB) 


GET REMAINING PARAMETERS 
(IACTI) 


READ AND PRINT ALLOY 
(BATCH I ) 


ACCESS ALLOY DATA BASE 
(ADB) 


READ AND PRINT REMAINING PARAMETERS 
(BATCH I ) 


EXTRACT PARAMETER VALUES 
FROM INPUT ARRAYS 
(INCON) 


CONVERT TO CGS UNITS 
(INCON) 


INPUT ERRORS 
IN THIS 
\ CASE?^ 


RETURN 











3. ) .*< Output Section 
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3.2 ALPHABETICAL LIST Of SUBROUTINES 

The list below contains only the modules that were written specifically for the 
solidification model. Off-the-shelf routines used by the model are listed In 
Sections 3.2 and 3.3 of reference 1, The purpose of each routine Is described 
briefly below; details of the programmlng» Including calling sequence descriptions, 
are in the program comments. 



NAME 

SECTION 

FUNCTION 

1 

ADB 

I/O 

Retrieve phase diagram, densities and viscosity from 
alloy data base. 

i| 

AD ISP 

I/O 

Displays any array on the screen or In printed form. 

li. 

ALNU 1 N 

I/O 

Gets Interactive Input that may be either alphanumeric 
or numerf cal . 

11 

AXES 

Graphics 

Draws and labels axes on all plots. 

U 

BATCH 1 

I/O 

Controls batch mode card Input. 

y 

COORD 

Calculation 

Calculates Isotherm positions and coordinate 
transformation. 


CTR 

Graphics 

Controls contour plots. 

II 

II 

ii 

CTRBOX 

Graphics 

Logic for contouring over given rectangular region. 

CTRDRW 

Graph i cs 

Draws contour line elements. 

rt ■ 

1! 

|| 

CTRLAB 

Graphics 

Labels contour level values. 

u 

DISKC 

I/O 

Performs various control functions on disk files. 

rr~- 

n 

y ■ 

DISK! 

I/O 

Retrles^es final solid data from temporary disk for 
use by graphics routines. 

n'' 

DISKO 

I/O 

Writes final solid data on temporary disk. 

ii 

ii 

EVAL 

Calculation 

Evaluates specified functions at current time. 

GPBLK 

Graphl cs 

Puts parameter block on plot. 

1 f 

GPHCON 

Graphics 

Controls interactive graphics. 

U ■ 

GRIDCN 

Calculation 

Controls size of computational grid. 

r 

GRIDXE 

Calculation 

Sets up computational grid. 

y. 

HPROFS 

Graph 1 cs 

Controls horizontal profiie plots. 


lACTI 

I/O 

Controls interactive input. 


INCON 

I/O 

Controis input section. 


INIT 

Calculation 

Initiai ization routine. 

I 

MASSCK 

I/O 

Monitors conservation of mass. 


MSG 

I/O 

Puts any brief message on terminal screen. 








NAME 

SECTION 

FUNCTION 

PERM 

Calculation 

Ca 1 cu 1 ates permeabl 1 1 ty . 

PSETUP 

Calculation 

Calculates t, and B and sets up boundary 
conditions (see equations 1.3.^.^ to l.3*^*n) 

PSOLVE 

Calculation 

Solves the pressure equation (see reference 1, 
Vol. 1, i<.5) 

SCALE 

I/O 

Calculates neat scales for plots. 

SETHOL 

1/0 

Sets up long hollerlth arrays. 

SETPLT 

Output 

Initializes all plots and puts on labels. 

START 

Calculation 

Takes first step away from chill. (See Vol. 1, 
Section k. 3 ) 

STEP 

Calculation 

Advances and 1 one step (See Vol 1, Section k.h) 

TABCON 

I/O 

Controls tabular output 

TMCON 

I/O 

Controls display checkpoints and time-varying I/O. 

VECPLT 

Output 

Controls plots of vector fields. 

VLCTY 

Calculation 

Calculates velocity. 

VPROFS 

Output 

Controls vertical profile plots. 

WAIT 

Output 

Waits for the operator to enter a P. 

3.3 key 

PROGRAM SYMBOLS 


A description of all common block items is at the end of subroutine INIT. 

PROGRAM 

SYMBOL 

COMMON 

BLOCK 

DEFINITION 

(See Volume 1, Section 2) 

AE 

/PROCSS/ 

^E 

AGC 

/SLZONE/ 

1 

ALNAM 

/ALLOY/ 

Element 1 is the hollerlth solvent name. 
Element 2 is the hollerlth solute name. 

BE, BL 

/PROCSS/ 

bj, b^ 

CL 

/SLZONE/ 


CL0 

/ALLOY/ 

C 

0 

DCLDT 

/ALLOY/ 

dC|_/dT 

DRHDC 

/ALLOY/ 

dpL/dCb 

DTDTM 

/SLZONE/ 

3T/at = e 

DTM 

/MESH/ 

At 

DXI 

/MESH/ 

A? 

DET 

/MESH/ 

An 

EPR 

/ALLOY/ 

k 
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PROGRAM 

SYMBOL 

COMMON 

BLOCK 

DEFINITION 
(See Volume 1 

ET 

/MESH/ 

n 

FLAG 

/SLZONE/ 

CjCxj) 

GAMMA 

/PMBLTY/ 

Y 

GL 

/SLZONE/ 

9l 

GMY 

/PROCSS/ 

9 

HMOLD 

/PROCSS/ 

L 

KK 

/PMBLTY/ 

K 

Nl 

/MESH/ 


NJ 

/MESH/ 


P 

/SLZONE/ 

P 

QLE 

/PROCSS/ 

q 

RHL 

/SLZONE/ 

Pl 

RHLE 

/ALLOY/ 

Ple 

RHS 

/ALLOY/ 

Ps 

RHSE 

/ALLOY/ 

PSE 

T 

/SLZONE/ 

T 

TE 

/ALLOY/ 

Te 

TM 

/MESH/ 

t 

V 

/SLZONE/ 


WMOLD 

/PROCSS/ 

""c 

X 

/MESH/ 


Y 

/MECH/ 

y 

XE, XL 

/PROCSS/ 


XI 

/MESH/ 



3.^ PROGRAM CONFIGURATION ON THE PRIME 


The procedures described In this section are used to maintain the program on 
the Prime ^00 system. 

3 .A.I Compilation 

If the FORTRAN source code is stored in a file named MPS2.PGM, then it can be 
compiled by entering the command 


FTN MPS2.PGM 2/500 


3.^.2 Unking 

The Prime utility for Unking end running segmented programs In SEG. It can 
be used to build a run file as follows: 

SEG 

LO ^MPS2 

LO B_MPS2.PGM 

LIB VAPPLB 

LIB TCS500 

LIB 

SAV 

Q 


3.5 EXECUTION COMMAND FILES 


After the run file #MPS2 has been built, the program can be executed by entering 
the commands described In Section 2, The execution Is set up and controlled by 
two command files listed below: 


Command File MPS2 


Command File MPS2. BATCH 


OPEN I CARD 1 1 
CO -END 


OPEN CARDS 1 1 
SEG ^MPS2 
C 1 2 3 
CO -END 


CARDS Is a disk file containing the batch card Input, and ICARD contains the 
single word INTERACTIVE. CARDS or ICARD Is accessed by the program via FORTRAN 
logical unit number 5- Other files used by the program are Ml.D.B, the data 
base file; DISK, the temporary disk storage file; and PRINT the batch printed 
output file. Ml.D.B Is accessed via logical unit number 7, DISK Is on logical 
unit number 8, and PRINT is on logical unit 6. No printed output Is generated 
by an interactive mode run* 
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SECTION k 
SAMPLE CASES 


INTERACTIVE CASE 

The case shown here Is the Al-*«.!)| Cu case that the interactive model will run 
If none of the Input parameters are changed from their default values. This 
case demonstrates the use of the checkpoint control to run the case through 
completion, and several output options that are new with the unsteady model. 

Any of the non- vector functions can be displayed on a contour plot, the one 
shown is contours of final composition over the entire Ingot. Note that the 
location of the function maximum and minimum are marked on the plot by an x 
and an o. The base level, Cg «* k.5% Cu in this case, is drawn with a heavy 
line., levels below the base level are drawn with dashed lines and those above 
are drawn with solid lines. Every fifth level is drawn with a heavy line. The 
selection of the base level and the interval between levels is under user control 
as shown on the plot selection page preceeding the sample contour plot. The 
next plot is horizontal profiles of final composition, using the default scaling 
showing the selection of three profiles. The last plot Is vertical profiles 
of final composition. Illustrating the change In plot scale. In addition to 
the plots shown here, all plots that were available In the steady-state model 
are also available In the unsteady model. The final output for this case is 
the tabular display of the metal mass balance and the solute mass balance at 
the end of solidification. 


n^TERlALS PROCCSSINO IN SPACE 
WCROSEGRECftTIOH IN A CASTING INGOT 
N<DEL 8 

1 UNIDIRECTIONAL SOUDIPICATION Of A IINARV f ILOV 
I UNSTEADV SOLUTION 

1 PLANAR ISOTHERMS, RECTANGULAR MUSHY ZONE 
1 ISOTHERM MOVEMENT INPUT 
I TEMPERATURE GRADIENT ASSUMED 
1 NO CONUECTION IN lULK LIQUID 
I ISOTROPIC PERMEAIILITV GAMMAtGL»8 


EN1ER CASE TITLE (UP TO 8« CHARACTERS) 


ALLOY 

50LUfN.T« AL 

SOLUTCI CU UCICKT PERCE.ITl 4.580E 00 

ENTER A TO CHANCE ALLCV OR P TO PROCEED. 

P 


ALLOY DATA lASE - DATA BASE REUISEO •4/30^80 


SOURCE OF IHFORrATION FOR AL -CU 0.O0eE-Ol TO 3.3C0E 01 UT. PCT. CU 

PHASE DIAGRAM - fl.C. FLENIhCS AMD C. MEREO, NET TRANS. UOL 335. 1567. P H<9. 
DENSITIES - R. MEhRABIAM. n. KEANE AND N.C. FLENIhCS. TRANS TNS-AINE, 

UOL 1. 1970. P 1309. 

UISCOSITV - R. NEHRABIAN. II. KEANE AND N.C. FLENIhCS. TRANS T.TS-AINE. 

UOL 1. 1970. P 1309. 


PHASE DIAGRAM 

TEMPER ATURE-CO.-'POS IT ION SLOPE 
EQUILIBRIUM PARTITION RATIO 
EUTECTIC COMPOSITION 
EUTECTIC TEMPERATURE 


DENSITIES 

COMPOSITION-DENSITY SLOPE 
SOLID DENSITY 
LIQUID EUTECTIC DENSITY 
SOLID EUTECTIC DENSITY 


UISCOSITV 


-3.89OE-01 PCT SOLUTE / DEC C 
I.730E-01 

3.30OE ei PCT SOLUTE 
5.A80E 03 DEC C 


a.G70E-C3 <C.VCnxt3) / PCT SOLUTE 
2.630E 00 Cn/Cnti3 
3.330E 00 cn/cni*3 
3.430E CO CM/CM113 


3.e0CE-03 Cn/(CMiSEC) 


ENTER P TO PROCEED 

P 
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SOLIDIFlCftTION rtODEL 


•9II9III THU. JAN 89 1981 


SOLIDIFICATION PROCESS PARANETERS 
1 HOLD (HALF) WIDTH B.7IIE M (CN) 

B HOLD HEIGHT S.SSfE •• (CN) 

3 ISOTHERM TINE EXPONENT (Q) 1.M8E •• 

A LIOUIDUS ISOTHERM COEFFICIENT (ID 1.7I8E-IB (CM/SEC) 

6 EUTECTIC ISOTHERM COEFFICIENT (IE) 1.7IIE-II (CM/-8EC) 

6 INITIAL ISOTHERM SEPARATION (AE) -B.tHC If (CM) 

7 CRAUITATIONAL FORCE l.MIE It (0) 


ENTER ITEM NUMBER TO CHANCE. OR 
P TO PROCEED. 

P 


AL 4. SIM CU SOLIDIFICATION MODEL 8 I9il9ill THU. JAN 89 1981 


PERMEAIILITV MODEL PARAMETERS 
1 GAMMA 


S.MIE-17 (CMttB) 


ENTER ITEM miMIER TO CHANGE. OR 
P TO PROCEED. 


4.5NI 


SOLIDIFICATION MODEL B 


19119118 THU, JAN 89 1981 


NUMERICAL METHODS CONTROL PARAMETERS 
1 NUMBER OF HORIZONTAL MESH POINTS 
B NUMBER OF UERTICAL MESH POINTS 
3 MAXIMUM NUMBER OF PRESSURE ITERATIONS 
A PRESSURE CONUERGENCE CRITERION 

5 FRACTION OF CFL STEP TO TAKE 

6 MAXIMUM STEP SIZE IN XE 

7 INITIAL UOLU)« FRACTION LIQUID AT CHILL 

8 MINIMUM HCmiZONTAL MESH SIZE 


ENTER ITEM NUMBER TO CHANGE, OR 
P TO PROCEED. 

P 


at 

II 

ail 

l.lllE-14 

6.IIIE-I1 

I.IIIE-Il 

9.IIIE-I1 

A 






CAICUIATIOM IN PROORESS FOR CASE 


AL 4.5NI CU 


SOLIDIFICATION NODEL 8 


•OtlOlll THU. JAN 89 ISIl 



TINE STEP SIZE 

PRESSURE CVCU6 

XE 

XL 


i.8>B34 I.MINI 

11 

>1.8€tll 

•.31890IE-I1 

ENTER 

T TO DISPLAY TABULAR DATA, 
G TO DISPLAY GRAPHS, 

C TO CONTINUE THIS CASE, 





0 TO TERNINATE RUN, OR 




C 

N TO PROCEED TO NEXT CASE. 





TINE SINCE 8EGINNING OF SOLIDIFICATION - 1.8BB34 SEC. 

LIQUIDU5 ISOTHERN IS •.31B998E->M CN FRON CHILL. 

EUTECTIC ISOTiHERN IS -l.flSSM CH FRON CHILL. 

ENTER EITHER THE TINE OF THE NEXT CHECKPOINT IN THE FORNAT T'UALUE 

OR THE LIOUIDOS LOCATION AT THE NEXT CHECKPOINT IN THE FORNAT L-UALUE 
OR THE EUTECTIC LOCATION AT THE NEXT CHECKPOINT IN THE FORNAT E»UALUE 

E-a.7 


NEXT CHECKPOINT UIIL IE AT T • B76.47I SEC.. XL > 4.6D68I CN.. XE - 8.69801 CN. 


TINE 

STEP SIZE 

PRESSURE CYCLES 

XE 

XL 

8.82351 

0.941170 

10 

-1.958N 

0.479997E-01 

4.83S87 

1.41176 

17 

-1.9SIH 

I.71B996E-01 

5.89418 

1.05882 

15 

-1.61N0 

0.:.'98993E-01 

6.617S6 

1.32352 

14 

-1.81750 

1. 118499 

8.87810 

1.65440 

23 

-1.8593I 

•.140684 

9.B50eS 

1.37867 

22 

-1.83504 

0.164061 

11.2591 

1.60844 

21 

-1.10860 

•.191405 

13.1356 

1.87658 

85 

-1.77669 

0.883305 

14.7776 

1.64195 

% 

26 

-1.74178 

0.861811 

% 

E4s!sai 

2. 01039 

•i* 

28 

B.n'ses 

4.17306 

247.381 

1.85952 

28 

8.B9546 

4.B0647 

849.101 

1.719B9 

87 

8.83470 

4.83470 

850.958 

1.85784 

88 

B.B66B8 

4.86688 

858.653 

1.69482 

88 

8. 89610 

4.89610 

254.199 

1.54596 

84 

8.38138 

4.3B13B 

255.894 

1.69484 

85 

8.35019 

4.36019 

257. Ml 

1.51577 

86 

8.37G95 

4.37B96 

258.765 

1.35549 

81 

8.39901 

4,39900 

260.284 

1.51878 

88 

8.48488 

4.48488 

261.611 

1.31729 

82 

8.44781 

4.44781 

268.744 

1.14233 

14 

8.46663 

4.46663 

264.070 

1.32634 

IS 

8.48918 

4.4B918 

265.161 

1.09010 

16 

8.6I770 

4.60770 

266.056 

0.895791 

9 

8.6BB94 

4.6BB94 

267.169 

1.11380 

10 

8.54187 

4.64187 

267.982 

0.812863 

10 

8.56569 

4.66669 

268.575 

0.592910 

10 

8.56677 

4.66677 

269.008 

0.432339 

10 

8.67311 

4.6T311 


SOLIDIFICATION CONPLETE. 


ENTER T TO DISPLAY TABULAR DATA, 
G TO DISPLAY GRAPHS, 

Q TO TERNINATE RUN. OR 
H TO PROCEED TO NEXT CASE. 
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•KRUNAL PAGE IS 
'(•' P<M)R QUAIJTY 






AL 4.5MI CU SOLIDiriCATION MODEL I •BilSill THU, JMN IB IBIi. 


FUHCTIOH TO MIOT 


ittmttit 


S/L ZONE 


itiiiitm 


l UCLOCITV 
i PRCSSUREl f-H 
3 HtCSSUAC - BULK HYDMSTHT2C M 
1 FRACTION LIQUID 
5 NASS FLOU 
fi SOLUTE FLOU 

imtmti INOOT iiitttuit 

7 FINAL LOCAL AUCRAGE CONPOSITION 
I UOLUNE FRACTION EUTECTIC 


ENTER ITEM HUNIEH OF FUNCTION TO RLOT, OR 
P TO PROCEED* 

7 

ENTER ITEM NUHIER OF PLOT TYPE 
4 


PLOT TYPE 

1 UERTICAL PROFILES 
I HORIIONTAL PROFILES 

3 UlCTOR FIELD 

4 CONTOURS 


FUNCTION NlNINUn - Q.774 FUNCTION NAXIIWN • S.3t 

THERE UILL BE I LEUELS SPACED AT INTERVALS B.SBMM ABOUT A BASE LEVEL OF 4.fBRSB 

ENTER r TO PROCEED UITH TNESE LEVELS, OR 
C TO CHANCE LEVELS 
C 


ENTER BASE LEVEL 
4.S 

ENTER LEVEL INTERVAL 
*2S 


FUNCTION NININUH • B.774 FUNCTION NAXIfWN ■ S.3B 

THERE UILL BE 17 LEVELS SPACED AT INTERVALS B.25MM ABOUT A BASE UVCL OF 4. MSN 

ENTER P TO PROCEED UITH THESE LEVELS, OR 
C TO CHANCE LEVELS 


P 




AL A.tm cu 


soiKiric^TioN nocci i 


•BlIBIII THU, JAM IB IBIl 


FUNCTION TO FLOT 


imiiiiii 


S/L ZONC 


ttiiiiiiM 


I ucLocnr 
I FFCSSURCi F-N 

3 FRCSSUFC - IUC.K HVBiWfTATIC F 

4 FFFCTION LIQUID 

5 IMSS FLOU 

I SOLUTE FLOU 

ifitiiitti IH 60 T Mtimm 

7 FINAL LOCAL AUCFA6E CONFOSITIOH 
I UOLUNC FRACTION EUTECTIC 


FLOT TVFE 

1 UERTICAL FROFILEI 
I HORIIOMTAL FROFILES 
9 VECTOR FIELD 
4 CONTOURS 


ENTER ITER NUMBER OF FUNCTION TO FLOT, OR 
F TO FROCEED. 

7 

ENTER ITER NUMBER OF FLOT TVFE 
I 


RINIRUH FUNCTION VALUE IS R.7740M 
RANIRUN FUNCTION VALUE IS S.29I9I 

AUTOMATIC SCALING VIELOS FLOT RANGCl I.SM TO S.StR UlTH REMOTE EKFONENT 

OVER > MAJOR TIC INTERVALS 

ENTER P TO PROCEED WITH AUTOMATIC SCALING. OR 
LOUER BOUND OF FLOT INTERVAL. 

F 


ENTER V VALUE OF PROFILE OR F TO FROCEED. 


ENTER V VALUE OF PROFILE OR P TO FROCEED. 
3.17S 


ENTER V VALUE OF PROFILE OR P TO PROCEED. 
S.3S 


ENTER V VALUE OF PROFILE OR P TO PROCEED. 
P 


ME COIWOSITIQN 




k-s 



§.m UITH iEMTC EKNfICKT • 


• a.7«Mt CN. raofl chill. 



SOLIBIFICATION MBCL I NltBUI THU* J»H IB IBIi 
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VnVt I S TMMILM MT«* 


« TP ifiPuV aUNTMt* 

P TO ftPWIjW Tf PUN, OP 
P TO PPOCm TO POCT o 


ctm. 


PL 4 .IHI CU SOLlBinCPTlON POKL G OSiOBlII THU. JPN IB IBlI 


1 PASS IPLl^^CF. TABLES 
I FIHAL LOCAL M;EPAQE COPPOSITIOH 

3 FPACTION EUTECTIC 

4 UOLUAE FRACTION LIQUID 

E LOCAL AVERAGE SOLID COPPOSITIOH 

5 VELOCITY 

7 PRESSURE (P-PO) 

• PRESSURE - BULK HYDROSTATIC P 
B TEPPERATURE 
11 LIQUID COPPOSITIOH 
11 LIQUID DENSITY 
le PRESSURE EQUATION COEFFICIINTS 


ENTER ITEP (ftlPBER OF TABLE TO DISPLAY. OR 
P TO PROCEED. 




li 


tT 

ti 

11. 


/r~ 


If 

a 


f 

i 


.a 


k-]\ 



cu 


MLIlirXOMTIM nONL I 


•KMiM ICt. KC M II 




*1.2 BATCH CASE 


This example shows the same Al-^i.gl^Cu case as example except It Is for 
zero gravity and It Is run In batch mode. Note that output Is at checkpoints 
only and only of current function values. Thls,> to get a table of final 
composition across the ingot the user must set up a checkpoint card at each 
value of Xg where data is required. The deck used to run this case Is shown 
below. In this there are three checkpoints. At the first final composition 
and volume fraction solid across the S/L zone when “ 0 (chill face) are 
selected by putting one's In columns 1? and 19. At the second only final 
composition Is printed, and at the third final composition and fraction eutectic 
are printed. 


lATCH 

BATCH EXAMPLE 
AL 
CU 
4.5 
2.7 
fi.35 
1. 

.#17 

.•17 

- 2 . 

I, 

6.E-7 

21 

3 

2ee 

l.EH 

.5 

.1 

.5 

4 

E •. 

E .5 

E 1. 

END 


1 1 
1 

11 


HOLD HALF UIDTH (CHI 
HOLD HEIGHT (CH) 

Q 

BL 

BE 

AE 

CRAblTV <G) 

GAHHA (CnXtE) 

NUHBER OF HORIZONTAL HESH POINTS 
NUHBER OF UERTICAL HESH POINTS 
HAXIHUn NUMBER OF PRESSURE ITERATIONS 
PRESSURE COHUERGENCE BOUND 
FRACTION OF CFL STEP 
MAXIMUM STEP IN XE 
GLI 

MINIMUM NI 


4-13 




IHftTERlALS PROCESSING IN SPACE 
HACROSEGREGATIOH IN A CASTING INGOT 
MODEL a 

t UNIDIRECTIONAL SOLIDIFICATION OF A DINARV ALLOY 
t UNSTEADY SOLUTION 

I PLANAR ISOTHERMS, RECTANGULAR MUSHY ZONE 
t ISOTHERM MOUEHENT INPUT 
I TEMPERATURE GRADIENT ASSUMED 
t NO CONUECTION IN BULl; LIQU’!^ 

» ISOTROPIC PERMEABILITY K* GAMMA*GL*ta 

AL 4.5000 CU SOLIDIFICATION MODEL E 09i55«30 THU, JAN 89 1981 
BATCH EXAMPLE 

1 


ALLOY DATA BASE - DATA BASE REUISED 4/30/80 


SOURCE OF INFORMATION FOR AL -CU 0.000E-O1 TO 3.300E 01 UT. PCT. CU 

PHASE DIAGRAM - M.C. FLEMINGS AND G. NEREO, MET TRANS, UOL 839, 1967, P 1449. 
DENSITIES - R. MEHRABIAN, Ho KEANE AND M.C. FLEMINGS, TRANS TMS-AIME, 

UOL 1, 1970, 1E09. 

UlSCOSITY - R. MEHRADIAN, H. KEANE AND M.C. FLEMINGS, TRANS TMS-AIME, 

UOL 1, 1970, P 1809. 


PHASE DIAGRAM 

TEMPERATURE-COMPOSITION SLOPE 
EQUILIBRIUM PARTITION RATIO 
EUTECTIC COMPOSITION 
EUTECTIC TEMPERATURE 


*a.B8tiE-01 PCT SOLUTE / DEO C 

t , 72 Ar.Ai 

3.300E 0) PCT SOLUTE 
5.480E 08 DEG C 


DENSITIES 

COMPOSITION-DENSITY SLOPE 
SOLID DENSITY 
LIQUID EUTECTIC DENSITY 
SOLID EUTECTIC DENSITY 


8.670E-0a (GH/CHtl3) / PCT SOLUTE 
8.680E 00 GM/CMSX3 
3.880E 00 GM/CHtt3 
3. 4008 00 GM/CMSI3 


UlSCOSITY 


3.000E-0E GH/(CHfS£C) 
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